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Abstract 

The gradient flow in non-abelian gauge theories on M 4 is defined by a local diffusion equation 
that evolves the gauge field as a function of the flow time in a gauge-covariant manner. Sim- 
ilarly to the case of the Langevin equation, the correlation functions of the time-dependent 
field can be expanded in perturbation theory, the Feynman rules being those of a renor- 
malizable field theory on M 4 x [0, oo). For any matter multiplet and to all loop orders, we 
show that the correlation functions are finite, i.e. do not require additional renormalization, 
once the theory in four dimensions is renormalized in the usual way. The flow thus maps 
the gauge field to a one-parameter family of smooth renormalized fields. 



1. Introduction 



The physics described by non-abelian gauge theories can be studied in many ways. 
Depending on the context and the questions to be answered, the desired information 
is extracted from correlation functions of local fields, expectation values of Wilson 
loops or the Schrodinger functional, for example. 

As explained in refs. [1,2], the gradient flow provides further opportunities to probe 
these theories and allows some of their otherwise elusive properties to be understood 
(in lattice gauge theory, the gradient flow is referred to as the Wilson flow, while in 
the mathematical literature it is commonly known as the Yang-Mills gradient flow). 
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Evidently, physically meaningful probes must be safe of ultra-violet divergences or 
must be such that these can be canceled by a well-defined renormalization procedure. 

In the case of the gradient flow, there is some evidence that the gauge field gen- 
erated by the flow does not require renormalization [1,2], but a formal proof of the 
absence of ultra-violet divergences at positive flow time has not been given so far. 
The aim of the present paper is to fill this gap through an all-order analysis of the 
flow in perturbation theory. With respect to the closely related case of the renor- 
malization of the Langevin equation discussed by Zinn-Justin and Zwanziger [3,4], 
there are two important differences, one being the absence of the noise term in the 
flow equation and the other the fact that the initial distribution of the gauge field 
(which is given by the functional integral of the theory considered) is not ignored. 

The perturbative analysis presented in this paper applies to renormalizable gauge 
theories with any compact simple gauge group and any matter multiplet. However, 
in order to simplify the discussion as much as possible, only the pure SU(iV) gauge 
theory with dimensional regularization will be considered, the generalization to other 
cases being straightforward (see sect. 8). 

In the following three sections, the Feynman rules for the correlation functions of 
the gauge field generated by the gradient flow are derived and are shown to be those 
of a local field theory with an extra dimension (the flow time). The finiteness of 
the correlation functions at positive flow time can then be established using power- 
counting and the BRS symmetry (sects. 6,7), but for illustration the divergent parts 
of a set of one- loop diagrams are first worked out in sect. 5. 



2. Iterative solution of the flow equation 

In this section, we introduce the gradient flow and derive the Feynman rules for the 
associated flow- line diagrams. The pure SU(iV) gauge theory in D = 4 — 2e euclidean 
dimensions is considered and the fundamental gauge field A fl (x) is normalized so that 
its action at bare coupling g is given by 




(2.1) 



duA v 



dvAn + [Ap, A v \ 



(2.2) 



(see appendix A for unexplained notation). 
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2.1 Definition of the gradient flow 

The gradient flow evolves the gauge field as a function of a parameter t > that is 
referred to as the flow time. Starting from the fundamental gauge field, 

B,\ t=0 = A„ (2.3) 

the time-dependent field B^it^x) is determined by the differential equation 

dtBp = D v G ViX + a D^B„ (2.4) 

GW = dpB,, - d v B^ + !>„ = 9„ + [B„, ■ ]. (2.5) 

The name "gradient flow" derives from the fact that first term on the right of eq. (2.4) 
is proportional to the gradient of the gauge action along the flow. Note that neither 
the initial condition (2.3) nor the flow equation (2.4) involve the gauge coupling. 

The second term on the right of eq. (2.4) is included in order to damp the evolution 
of the gauge degrees of freedom of the field. As in the case of the Langevin equation 
[4] , some technicalities in the perturbative analysis of the flow can be avoided in this 
way without affecting the evolution of the gauge- invariant observables. The latter 
are in fact independent of the parameter ceo, since the solutions of eq. (2.4) obtained 
at different values of ao are related by a (time-dependent) gauge transformation [1]. 

2.2 Expansion in powers of the fundamental gauge field 

Equation (2.4) may be split into a linear and remainder part according to 

8 t B^ = d v d v B„ + (a - \)d,,d v B v + (2.6) 
R„ = 2[B„,d„B„} - [B u ,d^B u ] + (a - l)[B^d u B„] + [B„, [B„,B„]}. (2.7) 
The linearized equation can be solved using the heat kernel 



K t (z), u = / —{(S^-ptfJe-* +p^e~^ }, (2.8) 

J p P 

where 

H(0- 
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Taking the boundary condition (2.3) into account, the flow equation may then be 
cast in the integral form 

B„(t,x) = J & D y{K t (x-y)^A v (y) + J ds K t _ s (x - y)^R„(s, y)}. (2.10) 

From this representation the retarded character of the equation is evident and it is 
also quite clear that the sensitivity to the initial value of the field dies away as t 
increases, although only slowly so at small momenta (ao is assumed to be positive). 
When passing to momentum space, 

B^t,x)= [ e^B^t,p), (2.11) 

the integral equation (2.10) becomes 

B„(t,p) = k t {p)^A v {p) + f dsk t _ s (p)^R u (s,p). (2.12) 

J o 

It is helpful at this point to introduce the vertices X^ 2 '^ and X^ 3 ^ through 

3 

K^P) = Y.h I ■ ■ [ (^) D 5(p + qi + ... + q n ) 

n=2 n ' J <K J 1" 

x xW(p, qi ,. . .,q n C\-::ilBll(t, -qi)-.. -q n ) (2.13) 

and the requirement that they are totally symmetric in the momentum-index combi- 
nations (qi, ai), . . . , (q n , v n ,a n ). The solution of the integral equation in powers 
of the fundamental gauge field, 

B«{t,p)=K t {p)^Al{p) + \ f dsK t _ s (p)^ I (27r) D 5(p-q-r) 

JO Jq,r 

x X (2 >°\p, -q, -r)l b ; a K s {q) p5 K s {r) aT A b MK{r) + ..., (2.14) 

is then obtained through iteration, i.e. by recursively inserting the equation on the 
right of itself. 

The vertices X^ 2,0 ) and X^ are given explicitly in appendix B. Note that the 
momentum-index combination (p,/j,,a) plays a special role in eq. (2.13). In partic- 
ular, the vertices are symmetric only in their other arguments. 
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Fig. 1. The diagrams contributing to B^(t,p) are directed tree graphs with a single 
external line (a little square is drawn at the end of this line). Flow lines always start 
from a one-point vertex (circle with a cross) or a flow vertex (circle) and end at another 
flow vertex if the line is not the external one. Each flow vertex has one outgoing flow 
line that corresponds to the momentum-index combination (p,fj,,a) in eq. (2.13). 



2.3 Flow-line diagrams 

The terms contributing to the expansion (2.14) in powers of the fundamental gauge 
field can be graphically represented by Feynman diagrams (see fig. 1). There are two 
kinds of vertices in these diagrams, the flow vertices X^ 2,0 - 1 and X^ 3 ' ^ introduced in 
the previous subsection and the one-point vertex 

= Alip). (2.15) 

P 

Each flow vertex is inserted at some flow time that is integrated from zero to infinity, 
while the one-point vertices reside at time zero. The vertices are connected through 
the flow lines 

t, \i,a s,v,b 

-i = 5 ab e(t-s)K t _ s (p)^, (2.16) 

P 

where t and s are the flow times at the endpoints of the line. In view of the retarded 
nature of the propagator (2.16), the flow time increases from zero at the one-point 
vertices to time t at the outer end of the external line as one follows the arrows 
in the diagram. In particular, the times associated with the vertices are effectively 
integrated only up to t (rather than infinity). 

As far the momenta, the index contractions and the symmetry factors are con- 
cerned, the Feynman rules are the usual ones. It is then not difficult to show that the 
sum of all diagrams solves the flow equation (2.12) order by order in the fundamental 
gauge field. 
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3. Perturbation theory 



The n-point correlation functions of the field B®(t, x) can be computed in perturba- 
tion theory by expanding the field in powers of the fundamental field, as explained 
sect. 2, and by calculating the correlation functions of the latter in the SU(7V) gauge 
theory as usual. We now show that the correlation functions can be directly obtained 
from a set of Feynman rules in D + 1 dimensions. 

3. 1 Gauge fixing 

The flow equation (2.4) is invariant under the infinitesimal transformation 

SB^ = D^co, (3.1) 

provided the time-dependence of w(t, x) G su(N) is such that 

d t uj = aoD^d^u. (3.2) 

Since the initial value of oj is unconstrained, these transformations generate the full 
gauge group at flow time t = and thus extend the gauge symmetry of the SU(iV) 
gauge theory to all flow times. 

The symmetry can be fixed as usual by including the gauge-fixing term 

% = / d D xtT{d^(x)d u A u (x)} (3.3) 
#0 J 

and the associated ghost action 

S c -c = -\ [ ^xtxid^D^x)} (3.4) 

in the total action of the theory, where c and c are the Faddeev-Popov ghost fields. 
As far as the n-point correlation functions of the fundamental gauge field and the 
ghosts are concerned, the Feynman rules are then the standard ones. Evidently, since 
the transformation (3.1) is an infinitesimal gauge variation, the expectation values 
of gauge-invariant expressions in the field generated by the flow are independent of 
the gauge parameter Aq. 
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3. 2 Gauge- field propagator 

The way in which the flow-line diagrams combine with the Feynman diagrams of 
the underlying theory is best explained by considering the two-point function of the 
time-dependent gauge field. To leading order in the gauge coupling, the flow-line 
diagram 

t, \i,a s,y ,b 

— «-<g> (3.5) 

p q 

is the only one that contributes to the correlation function. The contraction of the 
gauge fields at the one-point vertices then shows that 

(B;(t,p)Bl(s,q)) = (27T) D 5(p + q)5 ab g 2 D t+s (p)^ + O(g 4 ), (3.6) 

^(PW = p)2 {(W 2 - P^>~ tp2 + Aq W e " antp2 }- (3-7) 
This formula includes the mixed propagator 

(Al(p)Bt(s,q)) = (fy{0,p)Bt{s,q)) (3.8) 

as well as the two-point function of the gauge field at flow time zero. 

Since all three propagators are given by the same analytic expression, the same 
graphical symbol 

t, \i,a s,v,b 

s/WX/V = 5 ab g 2 D t+s (p)^ (3.9) 
P 

may be used for them. Note that the contraction of the one-point vertices always 
has the effect of converting the terminal flow lines to gauge-field lines. If one starts 
from the flow-line diagram 



t, \i,a 

' (3 - 10) 

P 

instead of the diagram (3.5), for example, the leading-order contribution is obtained 
by substituting the tree-level diagram for the correlation function of the gauge fields 
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at the three one-point vertices. This leads to a diagram 




(3-11) 



that has an ordinary three-point vertex (filled circle) and a flow vertex. Two of the 
lines attached to the latter are gauge-field lines instead of flow lines, but the expres- 
sion for the vertex is the same as before. 

3.3 Feynman rules in D + 1 dimensions 

The flow time will now be interpreted as an additional space-time coordinate. Since 
only non-negative times are considered, the D + 1 dimensional space is a half-space 
with a D dimensional boundary at flow time zero. The SU(iV) gauge theory lives 
at the boundary, while the field generated by the gradient flow extends to the extra 
dimension. 

From this point of view, the ordinary and the flow vertices represent boundary 
and bulk interaction terms, respectively, while the propagation of the fields in D + 1 
dimensions is described by the gauge-field propagator (3.9) and the flow propagator 
(2.16). The ghost-field propagator 

(c a (p)~c\ q )) = (2ir) D 5(p + q)6 ab g 2 D(p) + O( 5o 4 ), D(p) = 1 (3.12) 

on the other hand, is defined at flow time zero only. Graphically it is represented 
through a dotted line, 

" - b = 5 ab g 2 D(p), (3.13) 
P 

where the arrow is drawn in the direction from c to c. 

It should be quite clear at this point that the n-point functions of the field gener- 
ated by the gradient flow are given by Feynman diagrams in D + 1 dimensions with 
the graphical elements listed in the caption of fig. 2. Some special features of the 
Feynman rules are worth pointing out: 

(a) Each flow vertex has exactly one outward-directed flow line corresponding to the 
first momentum-index combination in eq. (2.13). The other lines attached to these 
vertices may be gauge- field lines or ingoing flow lines. 
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Fig. 2. Examples of diagrams contributing to the two- and three-point functions of 
the time-dependent gauge field. All diagrams are built from flow propagators (directed 
solid lines), gauge-field propagators (wiggly lines), ghost-field propagators (directed 
dotted lines) , flow vertices (open circles) and ordinary vertices at flow time zero (filled 
circles) . The little square at the end of the external lines indicates that they are not 
amputated. 

(b) Flow lines must start at a flow vertex and are either external or end at another 
flow vertex. The gauge-field lines, on the other hand, can start and end at both the 
flow vertices and the ordinary vertices. 

(c) The flow vertices are inserted at some flow time which is integrated from zero 
to infinity. All other vertices are at flow time zero. The flow times on which the 
propagators depend are the ones at the endpoints of the corresponding lines. 

(d) Diagrams with closed flow-line loops are set to zero. This rule derives from the 
fact that flow-line diagrams are tree diagrams and that the contraction of the gauge 
fields at the one-point vertices never leads to new flow lines. 



4. Field theory D+l dimensions 

The Feynman rules for the correlation functions of the time-dependent gauge field 
obtained in the previous section are those of a local field theory in D + 1 dimensions 
[4]. It is possible to show this through somewhat formal functional-integral mani- 
pulations, but one can also adopt a purely algebraic point of view, where the action 
in D + 1 dimensions merely serves as a generating function for the Feynman rules. 

4-1 Action 

Except for the boundary condition (2.3), the field will now be considered 

to be an independent field. One also needs to introduce a Lagrange-multiplier field 
Lfj,(t, x) = L®(t,x)T a with purely imaginary components. The action of the theory 
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in D + 1 dimensions is then given by [4] 



Stot = S+S g f + 5 C e + 5fl, (4.1) 
f'OO r 

S a = -2 J dtj d D xtr{Ly(t,x)(d t By - D v G V y - a Dyd v B u ){t,x)} . (4.2) 

In this framework, the flow equation (2.4) coincides with the field equation obtained 
by varying the action with respect to the Lagrange-multiplier field. Note that the 
latter is not required to satisfy any particular boundary conditions. 

4-2 Propagators 

The quadratic part of the action is the sum of 

J™ tit J d D xL*(t,x)(d t B« - d v d v Bl - (a - l)dyd u B^)(t,x) (4.3) 

and the quadratic part of the action in D dimensions. The bulk fields By, and Ly thus 
couple to each other only, but there is also an implicit coupling to the fundamental 
gauge field through the boundary condition (2.3). This complication can easily be 
overcome by substituting 

By(t,x) = J d D yK t (x-y)y U A u (y) + by(t,x). (4.4) 

The field by then satisfies homogenous boundary conditions, while the first term in 
eq. (4.4) solves the linearized flow equation and thus drops out in the action (4.3). 

It follows from these remarks that the propagators of the fundamental gauge field 
and the ghost fields coincide with the expressions given in the previous section. All 
other two-point functions vanish except for 

{bl (t, x)L b u (s, y)) | leadin§ order = 6« b H(t, x- s, y)^, (4.5) 

which is determined by the field equation 

{Sypd t - 5y P d a d a - (a - l)dyd p }H(t, x; s, y) pv = 8y V 8(t - s)S(x - y) (4.6) 

and the boundary condition Hy U (t, x; s, y)\ t=0 s>0 = 0. The unique solution of these 
equations is 

H(t, x; s, y)y V = 0(t - s)K t _ s (x - y)y V (4.7) 
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1 



2 



3 



Fig. 3. Examples of diagrams with flow-line loops. All such diagrams vanish as a 
consequence of the retarded nature of the flow propagator and the rules of dimensional 
regularization [4]. 

and the bL propagator is thus seen to coincide with the flow propagator. 

The two-point functions involving the B field may finally be calculated by combin- 
ing eq. (4.4) with the results obtained so far. Since the AL and the bb propagators 
vanish, one quickly finds that the BL propagator is equal to the flow propagator 
and that the BB and the AB propagators are equal to the expressions given in the 
previous section. 

4-3 Vertices 

The theory described by the action (4.1) has the vertices of the theory in D dimen- 
sions plus those deriving from the bulk action Sr. Recalling eqs. (2.6) and (2.13), it is 
straightforward to show that the LB 2 and LB 3 vertices generated by the interaction 
part of the latter, 



All vertices and propagators of the Feynman rules of the previous sections are thus 
recovered. A somewhat unusual feature of these rules is the presence of a field (the 
Lagrange multiplier) that propagates only through its mixing with the other fields. 
The algebraic structure of the Feynman rules is however entirely standard and does 
not involve any prescriptions apart from the ones deriving from the expansion of the 
action in powers of the fields. 

4-4 Flow lines and flow-line loops 

Flow lines represent the BL propagator. They can start and end at the flow vertices 
but are never attached to an ordinary vertex. Outward- and inward-directed external 
flow lines are external B and L lines, respectively. 




(4.8) 



coincide with the vertices X^- 2 ' ^ and X^ 3 '°\ 
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Diagrams with closed flow lines comply with these rules, but should be absent if a 
complete matching with the Feynman rules of the previous section is to be achieved 
(cf. point (d) at the end of sect. 3). Such diagrams are in fact equal to zero [4]. The 
diagrams 1 and 2 in fig. 3, for example, vanish because dimensional regularization 
sets the momentum integral to zero. If there are two or more vertices in the loop, as 
in diagram 3, the time integrations vanish, because the flow propagator is retarded 
and thus forces the flow times at the vertices to be squeezed to a range of measure 
zero (singularities in the time coordinates are excluded in view of the regularization 
of the momentum integral). 

If a lattice regularization is used, the vanishing of diagram 2 may not be guaran- 
teed, but one can always include a ghost field in the action that cancels the flow- line 
loops algebraically, i.e. at the level of the Feynman integrands [3]. With dimensional 
regularization, however, this device is not needed and is therefore omitted here. 



5. Sample calculation at one-loop order 



For illustration, the divergent parts of the one-loop diagrams that contribute to the 
gauge-field two-point functions are worked out in this section. The parameter and 
field renormalization required at flow time zero is then seen to cancel the singularities 
at all flow times. 

5.1 Computation of self- energy diagrams 

There are only few self-energy diagrams at one-loop order (see fig. 4) . The diagrams 
1 — 3 coincide with the ones contributing to the two-point function of the gauge field 
at flow time zero. We therefore merely quote the known result 

rS«| pole = <-<W -Wljp; (f - ^) (5.D 

for the sum of their singular parts. 

In the case of diagram 4, the Feynman integrand is a linear combination of terms 
of the form 

J ^^(Hri 2 } 2 e " t{ri2+#+Pf} ' U ' V G (5 ' 2) 
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Fig. 4. One-loop diagrams contributing to the AA (1—3), LA (4), LB (5, 6) and LL 
(7) vertex functions. The time-momentum-index combinations at the external legs on 
the left and right are (t,p,fjL,a) and (s,q,v,b), respectively (the flow times are absent 
in the case of the A legs). 

where k is the loop momentum and Q(k,p)^ u a homogeneous polynomial in k and 
p of degree 6 (the momentum-conservation factor (2n) D 5(p + q) is now always sup- 
pressed). This diagram is finite at all t > and in any dimension D, but eventually 
it will be multiplied by the external propagators and must then be integrated over 
t from to infinity. The behaviour of the integral small t thus matters and can give 
rise to singularities at D = 4. 

If f(t) is any smooth test function (an external propagator, for example), the time 
integral 

/"OO 

I(k,p)= / dt/(t)e-'^ fc2+t, ( fc+p ) 2 > (5.3) 
Jo 

can be worked out in an asymptotic series at large k, the leading term being 

The integration over the flow time t thus improves the degree of divergence of the 
momentum integral. In particular, the first term in the decomposition 

I(k, P ) = {l(k,p)--. ^ r^) + 7 ^ (5.5) 

v ' \ ' [u + v)(k + r) 2 J (u + v)(k + r) 2 y ' 
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(where r is an arbitrary external momentum) makes the integral convergent at D = 
4. One is then left with a logarithmic divergence and therefore a pole singularity 
equal to the one of the ordinary Feynman integral 



y >u + v J k (fc 2 ) 4 (fe + r) 2 ' 



(5.6) 



which is easily evaluated using standard techniques. 
Proceeding in this way, one obtains 

^(^Cl^^SrtWV^^^), (") 

rS(M,rt;t| pole = ^s m ,)s^ (-jL) , (5.8) 

for the divergent parts of the LA and Li? vertex functions. The diagram 7 and thus 
the LL vertex function are finite at D = 4. 

5.2 Renormalization 

The bare coupling and gauge-fixing parameter are related to the renormalized pa- 
rameters g and A through 

si = H 2e S f 2 Z, A = AZ 3 -\ (5.9) 

where fx is the normalization mass. To one-loop order, the renormalization constants 
Z and Z 3 are given by 

Z = l-^y+CV), i* = ^y, (5.10) 

* = ! + <V), C0 ^(H_i.), (M1) 

up to scheme-dependent finite terms. 

The renormalization of the fundamental gauge field, 

Al = Z^ZI /2 {A^ (5.12) 

14 



involves both renormalization constants as a result of the unconventional normaliza- 
tion conventions adopted in this paper. For reasons explained in sect. 7, the ghost 
fields are renormalized asymmetrically according to 

<? = ZsZ^Z^icnT, c a = Z^Z-^^r, (5.13) 

where Z 3 is the usual ghost renormalization constant. Note that eq. (5.13) is equiv- 
alent to the standard renormalization prescription for the correlation functions at 
flow time zero, because the fields c and c always occur in pairs and only the product 
of their renormalization factors matters. 

5.3 Do the bulk fields require renormalization? 

To one-loop order of perturbation theory, the question may be answered by explicitly 
calculating the divergent parts (if any) of bulk- field correlation functions. We first 
consider the two-point function of the B field and define its Lorentz-invariant parts 
A and B through 

(B;(t,p)Bt(s,q)) = (2ir) D 5(p + q)j-^ 

x {(S^P 2 -p„p v )A(t,s,p 2 ) + p ll p„B{t,s,p 2 )}. (5.14) 

All self-energy diagrams drawn in fig. 4 contribute to A and B. The diagrams 4 — 6 
actually make two contributions, because their external legs are different. 
In the renormalized perturbation expansion 

oo 

X = ^Y J 9 2l+2 X il \ X = AoiB, (5.15) 

1=0 

the leading-order coefficients are 

^(0) = e -(t+s) P 2 ? B (0) = A -l e - Q0 (i+ S ) P 2 _ (g J6) 

At the next order, the residues of the poles that derive from the renormalization of 
the coupling and the gauge-fixing parameter are thus given by 

res^ 1 )} =-6 ^ (0) , (5.17) 

Z factors 

res^ 1 )} = (co -b )B (0 \ (5.18) 

Z factors 
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Recalling the results obtained in subsect. 5.1, it is now straightforward to verify that 
these poles are canceled by the divergent parts of the one-loop diagrams. Up to this 
order of perturbation theory, the two-point function of the time-dependent gauge 
field is thus finite and does not require further renormalization. 

In the case of the BL, BA R , LA R and LL correlation functions, the cancellation 
of the singularities at D = 4 can be shown in the same way. There is actually little 
to prove in the last two instances, because these two-point functions vanish to all 
orders of perturbation theory (there are no diagrams with ingoing and no outgoing 
flow lines). All calculations reported in this section thus support the conjecture that 
the B and the L field do not need to be renormalized. 



6. BRS symmetry 

We now proceed with the general discussion that will lead to the proof of finiteness 
of the correlation functions of the bulk fields to all orders in the gauge coupling. As 
a first step, the BRS symmetry of the theory in D + 1 dimensions is reviewed in this 
section. 

The gauge fixing discussed in sect. 3.1 follows the standard procedure and there- 
fore leads to a theory with a BRS symmetry. This symmetry acts on the boundary 
fields in the usual way, but the transformation of the bulk fields requires the solu- 
tion of the diffusion equation (3.2) and consequently tends to be non-local. In their 
work on the renormalization of the Langevin equation, Zinn- Justin and Zwanziger 
[4] however showed that the locality of the transformation can be restored by intro- 
ducing additional ghost fields. 

6. 1 BRS transformation of the boundary fields 

The BRS variation of the unrenormalized fields A^, c and c is defined by [5,6] 

<fA„ = D^c, (6.1) 

5c = -c 2 , (6.2) 

5c = \ d^A ll . (6.3) 

Note that the components c a and c a of the ghost fields and the operator 5 anti- 
commute with one another. The product on the right of eq. (6.2), for example, is 
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given by 

c 2 = c a c b T a T b = l c a c bfabc T c^ (g^ 

and the Leibniz rule for the operator 5 must take its anti-commuting character into 
account. The action S + S g { + S C c and the measure in the functional integral of the 
theory in D dimensions are then easily shown to be BRS invariant. 

6.2 Bulk ghost fields 

The additional ghost fields d(t, x) and d(t, x) mentioned above live in D + 1 dimen- 
sions but are otherwise of the same kind as the Faddeev-Popov ghosts. Their action 
is 

S dd =-2j dtj d D xtr{d{t,x)(d t d-a D^d)(t,x)} (6.5) 
and the d field is required to satisfy the boundary condition 

d\t=o = c - ( 6 - 6 ) 

No boundary condition is imposed on the d field, which, in many respects, plays a 
role similar to the Lagrange- multiplier field in the case of the gauge field. 
The propagator of the ghost fields, 

(d a (t,p)~d b (s,q)) = (2n) D 5(p + q)5 ab 8(t- s)K t _ s (p) + O(g 2 ), (6.7) 

Ktip) = e~ a ° tp \ (6.8) 

can be worked out following the steps taken in subsect. 4.2. Graphically the propa- 
gator is represented by a dashed line, 

t,( l. . - S l b = S ab 6(t - s)K t (p), (6.9) 
P 

where the arrow is drawn in the direction from d to d. There is a mixed propagator, 
(d a (t,p)~c b (q)) = {2K) D 5{p + q)5 ab glb t (p) +O( 5o 4 ), (6.10) 

A(p) = 4 e-ao * 9 > ( 6 - n ) 
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as well, which coincides with the cc propagator at t = and is therefore represented 
by the same graphical symbol (a directed dotted line). The cd two-point function, 
on the other hand, vanishes to all orders. 
The action 

^interaction = " / dt f (^) D 5{p + q + r) 
JO Jp,q,r 

x (p, q, r); bc B^(t, -p)~d(t, -q) b d(t, -r) c , (6.12) 

also gives rise to a new flow vertex X^ 1,1 ) . In the Feynman diagrams, it is represented 
by an open circle as the other flow vertices. The explicit expression for the vertex 
is given in appendix B. 

As in the case of the flow-line loops discussed in subsect. 4.4, it is possible to show 
that diagrams with dd loops vanish. Note that ghost loops with mixed propagators 
do not exist, because there is no cd propagator and the vertices only couple c to c 
or d to d. The only non-zero ghost loops are thus the usual ones of the theory at 
flow time zero. 

6. 3 BRS variation of the fields in the bulk 

The BRS symmetry acts on the bulk fields according to [4] 



6B„ = D^d, (6.13) 

SL^lL^d], (6.14) 

Sd = -d 2 , (6.15) 

Sd = D^-{d,d}. (6.16) 



Note that 5B^ = SA^ and 5d = 5c at the boundary, as must be the case in view of 
the boundary conditions (2.3) and (6.6). 

In order to show that the BRS variation of the bulk action Sr + S d ^ vanishes, it 
is helpful to introduce the fields 

E„ = d t B„ - D v G vli - a D„d v B v , (6.17) 
e = d t d- a D^d^d. (6.18) 
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Fig. 5. Tree diagrams contributing to the correlation function on the left of eq. (6.23) 
(diagrams 1 and 2) and to the one on the right of the equation (diagrams 3,4 and 5). 
The momenta are ingoing and satisfy p + q + r — 0. In diagram 3, the momentum p 
flows into a vertex representing the insertion of the field [B^ , d] , while in the case of 
the diagrams 4 and 5 the external lines with momentum p are multiplied by ip^ (thus 
representing the field d^d). 

After some algebra, one finds that 



8E„ = [E^dj + D^e, 



(6.19) 



Se = — {e, d}, 



(6.20) 



and the invariance of the bulk action is then easily established. 
6.4 Examples of BRS identities 

Since the functional integration measure is also invariant, it follows that {SO) = 
for any product O of the basic fields. This leads to identities such as 



\ {d,Al(x)d v A b v (y)) = g 2 5 ab 5(x - y) 



(6.21) 
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Xo(B a Jt,x)d v A b v (y)) = -((D^dr(t,x)c b (y)), 



(6.22) 



\o(Bl(t,x)d v A b v (y)d p A° p (z)) = -((D fl dnt,x)c b (y)d p A c p (z)), 



(6.23) 



the first of them being the familiar gauge Ward identity in the theory in D dimen- 
sions, which determines the longitudinal part of the gauge-field two-point function. 

It is instructive to check eqs. (6.22) and (6.23) at tree level of perturbation theory. 
The first equation relates the longitudinal part of the mixed gauge-field propagator 
to the mixed ghost propagator, an identity that is immediate from the expressions 
for the propagators. The other equation involves the 3-point vertices (see fig. 5) and 
a non-trivial cancellation among various contributions (appendix C). 



7. Finiteness of the renormalized perturbation expansion 

The renormalization constants Z, Z% and Z3 are now assumed to be such that the 
singularities of the correlation functions of the fields {A R ) p ,c R and c R at D = 4 
cancel to all orders of the renormalized coupling. Our aim in this section is to show 
that the correlation functions involving the bulk fields at positive flow time are then 
finite too and thus do not require further renormalization. 

7. 1 Renormalized perturbation theory 

We first need to reorganize the perturbation expansion in terms of the renormalized 
parameters and fields. In this form, the expansion is generated by an action So + AS 
of the renormalized fields, which includes the counterterms AS required to cancel 
the poles of the Feynman diagrams at D = 4. 

Explicitly So and AS are obtained by expressing the bare parameters and fields 
in the total action 

Stot = $ + Sgf + Sec + S fl + S d j (7.1) 
through the renormalized ones and by setting 

So = Stot\ Z =Z 3 =Z 3 =l ' ( 7 - 2 ) 



AS = Si 



•tot — 



So + ASbc, 



(7.3) 
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where 

A5 bc = 2 J d D xtr{(Z^ 2 Zl /2 - 1)^(0, x)(A R ) M (x) 

+ {ZzZ^Zl' 2 - l)J(0,x)c R (x)}. (7.4) 
The boundary conditions are then 

B At=o = (^r)m' 4=o = C R> ( 7 - 5 ) 

and the Feynman rules derived from the action So thus coincide with those discussed 
in sect. 3 (apart from the fact that the bare parameters go and Ao are replaced by 
fi e g and A, respectively). 

Note that the boundary conditions (7.5) differ from the ones imposed on the bare 
fields (a product of renormalization factors is missing). The counterterm AS^c must 
be included in the action of the renormalized fields to correct for this. It amounts 
to adding two-point vertices to the Feynman rules, whose effect on the gauge-field 
and ghost propagators is equivalent to a change of the boundary conditions by the 
missing renormalization factors. 

7.2 Absence of bulk counterterms 

In the renormalized perturbation expansion, singularities at D = 4 (if any) appear 
for the first time at some loop order I > 1. Since the interactions are local, and 
since all propagators are tempered distributions in position space with singularities 
only at coinciding arguments, the divergent parts of the correlation functions at this 
loop order are expected to be such that they can be canceled by local counterterms. 
Similarly to the case of an ordinary field theory on a half-space studied by Symanzik 
[7], the counterterms can be localized either in the bulk or at the boundary of the 
half-space. 

In the theory considered here, bulk counterterms can be excluded from the outset. 
In order to show this, first note that the correlation functions of the bulk fields B^, 
L^, d and d are, at large flow times, given by Feynman diagrams built from flow 
lines and flow vertices only. All diagrams of this kind are directed trees or products 
such trees, each tree ending at one of the B and d fields in the correlation function 
considered. Moreover, the flow lines at the other ends of the trees must start from 
the L and d fields in the correlation function. 

Since there are no loop diagrams, the correlation functions of the bulk fields are 
non-singular at large flow times and do not require renormalization. Divergent bulk 
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counterterms are therefore excluded. Note, incidentally, that correlation functions 
of local fields composed from B and d fields are finite too, because no loop diagrams 
are generated when the arguments of some of these fields coalesce. The field D^d 
that appears in the BRS identity (6.23), for example, is of this kind and consequently 
does not need to be renormalized. 

7.3 Boundary counterterms 

The discussion in the previous subsection shows that the structure of any divergent 
parts of the correlation functions must correspond to the insertion of a counterterm 
localized at flow time zero. Since the correlation functions of the renormalized fields 
(^4r) m ,c r and c R are finite to all orders, divergences can only arise from diagrams 
with at least one flow vertex and thus at least one such vertex with an external flow 
line. The possible counterterms are therefore proportional to or d. 

Since the gauge coupling is dimensionless, the counterterm may not involve com- 
posite fields of dimension larger than 4 (such terms would be irrelevant). Moreover, 
Lorentz symmetry, global gauge invariance and the ghost number conservation must 
be respected. Since and d have dimension 3 and all other fields dimension 1, it 
follows that the possible counterterms at ^-loop order are of the form 

2g 21 J d D xtr{ Zl L p (0,x)(A R ) p (x) + z 2 d(0,x)c K (x)} (7.6) 

with some (singular) coefficients z x and z-i ■ An LB and a dd term should in principle 
be included here, but in view of the boundary conditions (7.5), these terms are not 
independent and their inclusion would be equivalent to a change of the coefficients 
z\ and Z2- 

7.4 Consequences of the BRS symmetry 

We now show that the boundary counterterm (7.6) is excluded by the BRS symmetry 
of the theory. In terms of the renormalized fields and parameters, the BRS identity 
(6.23) reads 

X(B;(t,x)d^A R t(y)d p (Anr p (z)) = -{{D,dY{t,x){c K ) b {y)d p {A K T p {z)). (7-7) 

The unusual asymmetric renormalization (5.13) of the ghost fields c and c was chosen 
to ensure that the renormalization factors drop out in this equation. From the point 
of view of the SU(./V) gauge theory, the standard and the asymmetric renormalization 
of the ghost fields are equivalent, but the situation is different in the theory in D + 1 
dimensions, because only c couples to the bulk fields. 
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Equation (7.7) holds as long as no extra counterterms need to be added and thus 
up to loop order I inclusive. Note that all fields in the correlation functions in this 
equation, including the composite field D^d, are renormalized fields that cannot be 
additionally renormalized (cf. subsect. 7.2). If the correlation functions are singular 
at Z-loop order, one must therefore be able cancel the singularities by a counterterm 
of the form (7.6). 

The contribution of the counterterm (7.6) to the correlation functions is obtained 
by inserting the corresponding two-point vertices in the tree diagrams in fig. 5. The 
insertions have the effect of multiplying the values of the diagrams 1 — 5 by g 21 times 

2z 1: zi, z 1 + z 2 , z\ + z 2 and z 2 (7.8) 

respectively. However, recalling the values of the diagrams quoted in appendix C, 
the sum of the diagrams weighted by the factors (7.8) turns out to violate the BRS 
identity (7.7) unless z\ = z 2 = 0. 

The correlation functions in eq. (7.7) must therefore be finite at Z-loop order and 
all other (renormalized) correlation functions must be non-singular too, because the 
addition to the action of a counterterm of the form (7.6) is excluded. We have thus 
shown that the theory in D + 1 dimensions does not require further renormalization. 



8. Miscellaneous remarks 



8. 1 Behaviour of the gauge field near flow time zero 

In the regularized theory, the time-dependent field B^{t, x) satisfies the boundary 
condition (2.3) and its correlation functions thus converge to those of the bare gauge 
field A^(x) when t goes to zero. However, since the latter requires renormalization 
by a divergent constant, the correlation functions tend to become singular at t = 
after renormalization and removal of the regularization. Their asymptotic behaviour 
for t — >■ is then described by an expansion 

B^x) = c B (t){A n )^x) + 0(t) (8.1) 

in local renormalized fields with singular coefficients [7] . At one-loop order of pertur- 
bation theory, for example, the coefficient ce(t) is found to diverge logarithmically 
at t = 0. 
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It is straightforward to show that the renormalization group equation 

r dji + %~ 2lX d\ + 7 + p/9 \ CB{t) = 0{t) (8 - 2) 

holds, where 

(3 = -b g 3 + O(g 5 ), 7 = C0 « 7 2 + O( 5 4 ), (8.3) 

are the beta function and the anomalous dimension of the gauge field. In the Landau 
gauge, for example, the equation can be easily integrated and one finds that 

B,(t,x) ~ o (2b g( q ) 2 ) 1/2 ~ Co/2bO R(9)(An),(x), (8.4) 

g(q) being the running coupling at momentum q = (8t) -1 / 2 and R(g) the factor that 
relates the renormalized to the renormalization- group-invariant gauge field. 

Equation (8.4) is a remnant of the boundary condition satisfied by the gauge field 
in the regularized theory. It shows that the field generated by the flow equation is 
connected to the fundamental field in a universal manner, i.e. there is no room for 
finite renormalizations here, the reason being that any such renormalization would 
violate the BRS symmetry. 

8.2 Gauge-invariant composite fields 

Wilson loops and gauge-invariant local fields are independent of the parameter i n 
the flow equation and do not require renormalization at positive flow time. At small 
flow times, their asymptotic behaviour is determined by the scaling properties of the 
corresponding renormalized fields in the SU(iV) gauge theory and thus reflects the 
singular nature of the latter. 

For illustration, consider the density 

E = jG^G^ (8.5) 

that has previously been studied in ref. [1]. Since E can mix with the unit field, the 
dominant term at small flow times, 

E(t,x) = {E(t,x))+c E (t){\F^F- v } R (x) + 0(t), (8.6) 

is its expectation value, while the first subleading term is proportional to the renor- 
malized action density of the fundamental field. The associated coefficients, (E(t, x)) 
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and ce(t), satisfy a renormalization group equation that allows their exact asymp- 
totic behaviour at small t to be worked out analytically. 

8.3 Theories with matter fields 

In the presence of matter fields, the gradient flow is defined by the same equations as 
in the pure gauge theory (subsect. 2.1). The interesting but rather more complicated 
case where the matter fields are included in the time evolution is not considered here. 

In a renormalizable theory, matter fields are scalar fields of dimension 1 or fermion 
fields of dimension 3/2. Provided the regularization preserves the gauge symmetry, 
our argumentation in sects. 6 and 7 then carries over literally. In particular, addi- 
tional boundary counterterms involving the matter fields are excluded by the global 
symmetries and the requirement that their dimension must be less than or equal to 
4. The finiteness of the correlation functions at positive flow times is therefore again 
guaranteed. 

QCD is a prominent example of a gauge theory that has all the required properties 
to ensure finiteness, but the same applies to many more theories of interest, including 
the SU(2) Higgs model, supersymmetric versions of QCD and technicolour theories. 

8.4 Lattice regularization 

While the Feynman rules are more complicated than with dimensional regularization, 
the perturbative analysis of the gradient flow on the lattice is not expected to run 
into fundamental difficulties. 

A possible choice of the gradient term in the lattice flow equation is the gradient of 
the Wilson action [1]. There is however no reason to choose this particular action or 
to require that it coincides with the gauge action of the theory. As long as the term 
has the correct form and normalization in the classical continuum limit (and thus at 
tree-level of perturbation theory), the correlation functions of the time-dependent 
gauge field will not depend on the exact choices one makes, except for lattice effects 
vanishing proportionally to a positive power of the lattice spacing a. 

The continuous Lorentz symmetry is broken on the lattice, but the remaining exact 
symmetries are sufficient to exclude counterterms that have not already appeared 
in the continuum theory. We therefore expect that the correlation functions of the 
time-dependent fields do not require renormalization and that their values in the 
continuum limit are independent of the regularization (up to finite renormalizations 
of the coupling and the gauge- fixing parameter). Moreover, if the lattice theory is 
0(a) improved, the correlation functions are automatically improved too, because 
there are no candidate 0(a) counterterms with all the required properties. 
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9. Conclusions 



The fact that the gauge field generated by the gradient flow does not require renor- 
malization is a consequence of the locality, the symmetries and the parabolic nature 
of the flow equation. In particular, in the associated field theory in D+l dimensions, 
bulk counterterms are excluded simply because the evolution of the bulk fields is 
retarded at large times and thus described by tree diagrams. The absence of bound- 
ary counterterms other than those needed for the renormalization of the theory at 
flow time zero is however non-trivial and can only be shown using power-counting 
and the BRS symmetry. 

In presence of matter fields, the situation is essentially unchanged as long as only 
the gauge field is evolved in time. The finiteness of the correlation functions of the 
field at positive flow time is therefore still guaranteed. Including all or some of the 
matter fields in the flow is however an interesting option that remains to be explored. 

We thank Jean Zinn-Justin for helpful discussions on the renormalization of the 
Langevin equation and for encouraging us to proceed directly with the flow equation 
rather than considering the latter to be a limit of the Langevin equation. 



Appendix A. Notational conventions 



The Lie algebra su(iV) of SU(iV) may be identified with the linear space of all anti- 
hermitian traceless N x N matrices. With respect to a basis T a , a = 1, . . . , N 2 — 1, 
of such matrices, the elements X € su(N) are given by X = X a T a with real compo- 
nents X a (repeated group indices are automatically summed over). The structure 
constants f abc in the commutator relation 

[T a ,T b ] = f abc T c (A.l) 

are real and totally anti-symmetric in the indices if the normalization condition 

tr{T a T b } = -\5 ab (A.2) 

is imposed. Moreover, f acd f bcd = N5 ab . 

Gauge fields in the continuum theory take values in the Lie algebra of the gauge 
group. Lorentz indices [i, v, . . . are automatically summed over when they occur in 
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matching pairs. The space-time metric is assumed to be euclidean. In particular, 
p 2 = for any momentum p and 8^ = D. 

Appendix B. Flow vertices 

The flow vertices X^°\ and XW are defined through eqs. (2.13) and (6.12), 

respectively. They are explicitly given by 

X^ \ P ,q,r)«% = if abc {(r - q)^5 up + 2q p S^ - 2rJ^ p 

+ {oLo-l)(q v 5 lip -r p 5 llv )}, (B.l) 

X^°\ P ,q,r,sTXi = f abe r de {5^5 vp - <5 W M 

+f ade f bce (S^ - 6^6 pa ) + f ace f dbe (S^S pa - 5^S up ), (B.2) 

X^ 1 \p,q,r)l bc = a Q if abc r tl . (B.3) 

In the Feynman diagrams, the momenta in these formulae are identified with the 
ingoing line momenta. 



Appendix C. Evaluation of the diagrams in fig. 5 

The contributions of the diagrams in fig. 5 to the Fourier transform of the correlation 
functions in the BRS identity (6.23) are of the form 

„4„- fabc 

where n is the diagram number. Explicitly, one obtains 

C,x,i = \p\q ~ r) M e- to0 ^ 2+r2 ) + \{q 2 - r 2 )p^- ta ^ - C„, 2) (C.2) 

C,,2 = \{p 2 (q ~ r), + (q 2 - r 2 )p,}^\ (C.3) 
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for the diagrams contributing to the left-hand side of the equation and 




? r e -t ao (i 2 +r 2 ) 



(C.4) 





} 



(C.5) 




(qr + r 2 )p^e 



(C.6) 



for the diagrams on the right-hand side. The BRS symmetry requires 



(C.7) 



which is indeed the case. 
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